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We study analytically, via the Newman-Penrose formalism, the late-time 
decay of linear electromagnetic and gravitational perturbations along the 
■ event horizon (EH) of black holes. We first analyze in detail the case of a 

. Schwarzschild black hole. Using a straightforward local analysis near the EH, 

we show that, generically, the "ingoing" (s > 0) component of the perturbing 
| field dies off along the EH more rapidly than its "outgoing" (s < 0) coun- 

terpart. Thus, while along r = const > 2M lines both components of the 
^ ■ perturbation admit the well-known t~ 2l ~ 3 decay rate, one finds that along 

the EH the s < component dies off in advanced-time v as v~ 2l ~ 3 , whereas 
the s > component dies off as t> -2 ' -4 . We then describe the extension of this 
Q\ \ analysis to a Kerr black hole. We conclude that for axially symmetric modes 

the situation is analogous to the Schwarzschild case. However, for non-axially 
symmetric modes both s > and s < fields decay at the same rate (unlike 
i ■ in the Schwarzschild case). 
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I. INTRODUCTION 



When a gravitational collapse results in the formation of a black hole (BH), the grav- 
itational field outside the the event horizon (EH) relaxes at late time to the stationary 
Kerr-Newman geometry. Also, when the (pure) Kerr-Newman field external to a BH is 
perturbed by gravitational or electromagnetic waves, the perturbing field die off at late 
time everywhere outside the BH, and along its EH. In both scenarios, it is implied by the 
"no hair 11 principle that when the BH geometry settles down into its stationary state, all 
characteristics of the initial state (or initial perturbation) must somehow be lost, except 
for the conserved quntities associated with it: its total mass, electric charge, and angular 
momentum. (For a detailed review of the "no hair" theorems by Hawking, Israel, Carter, 
and Robinson, see [|TJ.) 

Remarkable as the "no-hair" principle is, it still gives no information about the mech- 
anism through which this "compulsory" relaxation process occurs. For example, it tells us 
nothing about the rate of the decay process. Clearly, such a detailed description of the 
late time decay is important not only for gaining more insight into the "no hair" principle, 
but, more practically, by virtue of the recent prospects of detecting gravitational radiation 
from astrophysical black hole systems. Also, the characteristics of the decay along the event 
horizon has an impact on the nature of the singularity along the inner horizon of charged 
@ and rotating || black holes. 

A detailed description of the late time decay outside a Schwarzschild black holes was 
first given by Price (for scalar and metric perturbations |4[], and for all integer-spin fields 
in the Newman- Penrose formalism ||). Price found that any radiative multipole mode l,m 
of an initially-compact linear perturbation dies off at late time as t~ 2l ~ 3 (where t is the 
Schwarzschild time coordinate). If a static multipole mode existed prior to the formation 
of the BH, then this mode will decay as t~ 2l ~ 2 . Price found these power law decay tails to 
be the same for all kinds of perturbations, whether scalar, electromagnetic or gravitational 
(and in this respect, the scalar field model proved as a useful toy-model for the more realistic 
fields). 

Price's results were later confirmed using several different approaches, both analytic and 
numerical P-pf , and where generalized to other spherically symmetric spacetimes |7|JT2]-[T7| . 
The validity of the perturbative (linear) approach was supported by numerical analyses of 
the fully nonlinear dynamics |T8| , |T4|| , indicating virtually the same power law indices for the 
late-time decay. 

Recently, several authors addressed the issue of the late-time decay of fields outside 
rotating black holes. First, a numerical simulation of the evolution of linear scalar |1J and 
gravitational |2l| waves on the background of a Kerr black hole was carried out by Krivan 
et al. Later, an analytic treatment of this problem (in the time domain) was presented by 
Barack and Ori [f2~T|-[2"5(] (following a preliminary analysis by Ori |23J|). Then, a study of 
the late time decay in Kerr using a frequency- domain approach has been carried out by 
Hod, both for a scalar field P3 and for nonzero-spin Newman-Penrose fields |ZH] (following 
preliminary considerations by Andersson [[J). 

The above analyses all indicate that power law tails characterize the decay in the Kerr 
background as well. In this case, however, the lack of spherical symmetry causes coupling 
between various multipoles. As a result of this coupling, the power-law indices of specific 



2 



spherical-harmonics multipoles are found to be different, in general, from the ones obtained 
in spherically symmetric black holes. Another phenomenon caused by rotation (first observed 
in |[24||) is the oscillatory nature of the late time tails along the null generators of the EH of 



the Kerr BH for nonaxially-symmetric perturbation modes. (See [p2| , p3| for details.) 

As we just mentioned, power law tails are observed not only at timelike infinity, but also 
at future null infinity and along the (future) event horizon. Several authors have analyzed 
the late time behavior of a scalar field along the EH of a Schwarzschild BH [|7|,|rTJ and a 



Kerr BH |p2| , [25| . In both cases, the power law indices of the late-time decay along the EH 
were found to be the same as along any fixed-r world-line outside the BH (apart from the 
above mentioned oscillations along the EH in the Kerr case). Thus, in the Schwarzschild 
case, an l,m scalar perturbation mode is found to decay along the EH as v~ 2l ~ 3 (or v~ 2l ~ 2 
for an initially static mode). (Here, v is an advanced-time coordinate, which we define in 
the sequel.) 

Quite surprisingly, a careful and thorough study of the behavior of realistic physical fields 
(electromagnetic and gravitational) along the EH has not been carried out so far (to the best 



of our knowledge), even in the Schwarzschild case ||27||. One would expect the scalar-field 



model to provide, again, a reliable picture of the actual behavior of realistic physical fields; 
however, a careful analysis of the behavior of such realistic fields at the EH reveals several 
interesting new features, uncovered by the scalar-field case. These features arise already in 
the Schwarzschild case, and thus we find it instructive to study and explain this simpler case 
first. Accordingly, in this paper we first explore in detail the behavior of electromagnetic 
and gravitational perturbations at the EH of a Schwarzschild BH. Then we describe the 
extension of this analysis to the Kerr case, and derive the power-law indices at the EH. Full 
detail of the analysis of the Kerr case will be given in a forthcoming paper [23] (as part of 
a comprehensive analysis of the late time decay of perturbations in the Kerr spacetime). 

We shall apply a linear perturbation analysis, based on the Newman-Penrose formalism. 
In this framework, a single master equation governs the (gauge- invariant) radiative parts of 
the linear perturbations of both the Maxwell tensor and the Weyl tensor. For both fields, 
our analysis reveals that the "ingoing" (s > 0) part of the perturbing field dies off at late 
time along the EH of the Schwarzschild BH faster than its "outgoing" (s < 0) counterpart: 
Whereas the s < fields admit the usual v~ 2l ~ 3 law, the s > fields decay at the EH like 
v - 2l ~ 4 . In the Kerr case, the above difference in the behavior of the s > and s < fields 
occurs only for axially symmetric (m = 0) modes; for non-axially symmetric modes, one 
finds the same decay rates for both s > and s < 0. These results are summarized in Eqs. 
(p5|), (|96|), and ( |97D below, in the concluding section. We also comment there about the 
significance of our results to the study of the interior of spinning black holes, and discuss 
the relation of our analysis to previous works [ f27|j . 

An important role in our analysis is played by the static solutions of the field equation. 
These turn out to show a peculiarity: As in the scalar field case, there is a static solution 
regular at the horizon, and a second, independent, solution which is irregular there. However, 
for s > fields, regularity of a static solution cannot be judged merely from its leading- 
order behavior at the EH. Rather, the distinction between the regular and irregular solutions 
involves the identification of a certain, sub-dominant, logarithmic term in the latter. Another 
peculiarity has to do with the relation between static solutions and monochromatic solutions 
(i.e. modes of a single Fourier frequency uj). For s > fields, unlike the scalar field (and 
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unlike the s < case), the EH-regular static solution cannot be approached from an EH- 
regular monochromatic solution by naively taking the limit u — > 0. One finds that for s > 
this naive limit leads to a static solution irregular at the EH. We study these unexpected 
features, and then qualitatively explain them using a simple (scalar-field based) toy model. 

The paper is arranged as follows: In Sec. we give some definitions and notations, 
and briefly review the Newman-Penrose formalism for perturbations of the Schwarzschild 
geometry. In Sec. |TTT| we introduce the late time expansion, to be employed in our analysis. 
The static solutions to the field equation, central to our analysis, are obtained in Sec. [TV], 
followed (in Sec. [V]) by a formulation of regularity criteria for physical fields at the EH. This 
puts us in position to analyze (in Sec. [VTD the late time behavior of physical fields along the 
EH. This analysis yields the power index for both s < and s > fields. Another perspective 
on the subject is obtained in Sec. [VTTL where we consider the behavior of momochromatic 



modes. In Sec. [VIII| we then introduce a simple toy-model, which yields further insight into 
our results. The extension of our analysis to the case of a Kerr BH is described in Sec. 
1X| . In the concluding section (sec. jX|) we summarize the results and discuss their physical 



significance and their relation to other works. 

II. DEFINITIONS AND NOTATIONS 

The line element in the Schwarzschild spacetime reads, in the standard Schwarzschild 
coordinates t, r, 9, (f, 

ds 2 = -(A/r 2 ) alt 2 + (r 2 /A) dr 2 + r 2 (d6 2 + sin 2 6 dip 2 ), (1) 

where M is the mass of the BH, and 

A(r) = r 2 - 2Mr (2) 

is a function which vanishes at the EH, r = 2M. Here, and throughout this paper, we use 
relativistic units, c = G = 1. 

As this paper concerns with the behavior near the event horizon (EH), we shall find it 
convenient in the sequel to introduce a new (dimensionless) radial coordinate, 

r-2M 

z = 3 

2M 1 ; 

which vanishes at the horizon. 

We shall also need the EH-regular (Kruskal) null coordinates 

where v = t + r* and u =t — are the Eddington-Finkelstein null coordinates, with 

n = r + 2M\nz. (5) 

To discuss perturbations of the Schwarzschild BH via the Newman-Penrose formalism, we 
introduce the tetrad basis of null vectors (Z M , n^, m M , m* M ), defined in the (t, r, 6, ip) coordinate 
system by |]5|j29[ 
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F= [r 2 /A, 1,0,0] 
= [l,-A/r 2 ,0,0] /2 

= [0,0,1, if sin 6}/(2 1/2 r). (6) 

(The components of the forth tetrad leg, m* M , are obtained from the components of m M by 
complex conjugation.) 



In the framework of the Newman-Penrose formalism [30| the gravitational field in vacuum 
is completely described by five complex scalars, \I>o • • • ^4> constructed from the Weyl tensor 
by projecting it on the above tetrad basis. Likewise, the electromagnetic field is completely 
described by the three complex scalars ip ,ip 1 ,ip 2 , constructed by similarly projecting the 
Maxwell tensor. In particular, 

= -C^jfrn^Fm*, ^ 4 = -C^s^m^n^m* 6 (7) 

represent the ingoing and outgoing radiative parts, respectively, of the Weyl tensor, and 

^o = yV, ^ 2 = i^m^n" (8) 

represent the ingoing and outgoing radiative parts of the electromagnetic field. 

In the Schwarzschild (unperturbed) background all Weyl scalars but ^2 vanish (as di- 
rectly implied by the Goldberg-Sachs theorem, in view of the Schwarzschild spacetime being 
of Petrov type D; see Sec. 9b, c in ||31|| ). In the framework of a linear perturbation analysis, 
the symbols ^1, 8^/2, ^3, ^4 and <p , Sipi, ip 2 are thus used to represent first-order pertur- 
bations of the corresponding fields (with 8^2 = ^2 — ^ background^ e ^ c y Q ne can g^Q^y ( see 
Sec. 29b in pifl ) that \l/o and \l/4, and also (po and <p<z, are invariant under gauge transfor- 
mations (namely, under infinitesimal rotations of the null basis and infinitesimal coordinate 
transformations). The scalars ^1 and ^3 are not gauge invariant, and may be nullified by 
a suitable rotation of the null frame. The entities 8^2 and 8(fi represent perturbations of 
the "Coulomb-like", non-radiative, part of the fields (in fact, one can also nullify 8^2 by a 
suitable infinitesimal coordinate transformation.) It is therefore only the scalars defined in 
Eqs. (0) and © which carry significant information about the radiative part of the fields. 
(Note, however, that gauge invariance of the radiative fields is guaranteed only within the 
framework of linear perturbation theory.) 

There is a single master equation governing linear perturbations of both the gravitational 
and the electromagnetic radiative fields defined in Eqs. (|7|||) |32|. In vacuum, this master 
perturbation equation reads 

r*A-^ tt - A-(A*"*'), r - J-^mn^ - - 

sin 9 sin 9 sin 9 

-2s[Mr 2 /A - r]*f ( + (s 2 cot 2 9 - s)^ s = 0, (9) 

where ^ s (t, r, 9, (p) represents the various radiative fields according to the following substi- 
tutions: 

= ^ s =+ 2 

# 4 = r -^ s= - 2 . (10) 
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In Eq. (|9|), the angular dependence of is separable through a decomposition in spin- 
weighted spherical harmonics | 33fl . 



oo I 



l=\s\ m=—l 



(11) 



The time-radial functions ip slm (r,t) then satisfy the field equation 



rfy« m - A- s+i (A s+ V s D,r + 2sr 2 (r - 3M)^f + (Z - s)(Z + s + l)A^< m = 0. (12) 



—•5+1 / A s+1 ./.slm 



slm 



■slm 



III. THE LATE-TIME EXPANSION 

In order to analyze the power-law decay of perturbations at late time, we decompose the 
field in the form 



r lm (r,t) = j2 F k lm ( r ) 



-kg — k 



(13) 



fc=0 



to which we refer as the late-time expansion |34[]. Substitution in Eq. (|T!|) yields an ordinary 
equation for each function F, 



slm. 



where D is a differential operator defined by 

D sl = Ad 2 /dr 2 + 2{s + l)(r - M)d/dr — (I — s)(l + s + 1), 
and the source term Sf} m is given by 



S s k im ee 2(k + k-l)r 



d(rF£ l ™ 
dr 



+ 2srMA- 1 F° l ™ 



(14) 



(15) 



(16) 



(We take F£ l ™ = 0.) 

The dominant late-time decay at world-lines of fixed r is described by the term k = in 



Eq. (13). To the leading order in 1/v, we have 

^ slm (r,t) = F° l ™(r)v- ko . 

Substituting v = t + r*, we also find, to the leading order in l/t, 

^ slm (r,t) = F fc t™(r) t~ k °, 



(17) 



(18) 



which, using the well-known result by Price ||, implies k Q = 21 + 3 (or k = 21 + 2 if a static 
mode Z is initially present). [] 



For brevity, we hereafter consider modes without initial static multipoles. 
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Since Sf m vanishes for k = 0, the term F£ l J£(r) satisfies the homogeneous differential 
equation 

D sl (Fg lm ) = 0. (19) 

This is just the field equation of a static mode /, m. Thus, F^(r) must be a static solution 
of the field equation. In the next section we shall study the static solutions for focusing 
attention on their asymptotic behavior at the EH. 



IV. STATIC SOLUTIONS 

Since Eq. (|19|) is a second-order differential equation, its general solution is spanned by 
two basis solutions. We shall primarily be interested in the asymptotic behavior of these 
solutions near the EH. The leading-order asymptotic behavior can be easily obtained from 
the asymptotic form of Eq. (|T9|) near the EH: One finds that for both s > and s < 0, the 
two asymptotic solutions behave there like 

^m^ A o and ^slm^ A -s ( 20 ) 

(to the leading order in A). However, it is also possible to write an exact, global, basis of 
solutions to the static equation, as we do now. 

In terms of the radial variable z, Eq. (|T9"| ) takes the form 

(-z)(l + z)Fg + [-(a + 1) - 2(s + l)z] F^ + (l- s)(l + s + 1)F Q = 0, (21) 

where a prime denotes d/dz. This is the hypergeometric equation |35[] for F (—z). One 
solution for Eq. fl2~T|) is given by (see Sec. 2.1.1 in f35fl ) 



, ,s _ f F(-l + s,l + s + l;s + l;-z) = 0+ (for s > 0), 
<pr ^ j ~ \ {4M 2 z)- s F{-l, I + 1; -s + 1; -z) = <p; (for s < 0), 1 > 

where F denotes the hypergeometric function. (Hereafter we often omit the indices slm for 
brevity). Note that since in both cases the first index is a non-positive integer, F is simply 
a polynomial in z, and so is <p r . (We choose this notation because, as we shall see below, 
<f) r is physically regular at the EH, whereas the other static solution, to which we shall later 
refer as 0j r , is irregular there.) The normalization in Eq. ( p2|) was chosen so as to conform 
with Eq. ([2DD [recall that at z = 0, the hypergeometric function F = 1, and note also the 
relation A = AM 2 z(z + 1)]. Thus, to the leading order in A, <p T is given by 

A° (for a > 0), , . 

A- s (for s < 0). { } 



A second, independent, static solution is given by (see Sec. 2.2.2, case 21, in |35 



, (4M^)-(l + z)-'- 1 F[/- S + l,/ + l;2/ + 2;(l + ^)- 1 ]=0+ (for S > 0), 
' " 1 " ' 1 (1 - zy^Fll + s + 1, 1 + 1; 21 + 2; (1 + ^)- x ] = 07 (for s < 0), 

(24) 
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where A s i is a normalization factor, 



A 8l = 1/F(l -\ a \+l,l + 1; 21 + 2; 1) = (2f "^^l 1} , (25) 



(cf. Eq. (46) in Sec. 2.8 of [p5|), chosen such that 0j r takes the simple leading-order asymp- 
totic form at the EH, namely 

A careful study of the asymptotic behavior of 0j r at the EH reveals that it includes a 
(sub-dominant) logarithmic term.| To analyze this logarithmic term, it is instructive to 
express the irregular solution in terms of <p r via the Wronskian method. The Wronskian 
associated with the homogeneous equation ( fL9|) is 

W = A" 5-1 , (27) 

and thus a static solution independent of (fi r may be expressed as 

4> ir = -0 r (r) J <pf{r')W{r')dr' 

= -0 r (r) f <j ) ; 2 {r , )[^{r , )]- s - 1 dr' . (28) 



This solution, of course, does not necessarily coincide with ir , but is, in general, a linear 
combination of the two basis functions <p r and 4>i r . It is easy to verify that the integrand in 
Eq. (|28| ) is times a rational function which is regular (and nonvanishing) at z = 0. 

The integrand can therefore be expanded as 

z~ |shl (7o + liz + l2 z 2 + •••), (29) 

where 7» are constants, with 70 = (AM 2 )^^ 1 7^ 0. By substituting the leading-order term 
of this expansion in Eq. fl28|) and comparing to Eq. (|26|), we find 

ir = (2M|s|)"V<r + const • <f> r . (30) 

[Here, the coefficient of r depends on the specific choice of the lower integration limit in 
Eq. (p8[)-1 Substitution of the full expansion (p9|) in Eq. (|28|) yields 

ir = -<j> r [7| S | In z + ^" |s| (7o + 7i^H )] , (31) 

where 7, = 2M r y i /(i — \s\) for i 7^ |s| (and 7U1 is an arbitrary integration constant). We now 
use Eq. ( |30"D to extract ir : 



2 Such a logarithmic term is to be anticipated, because of the integer difference, |s|, between the 
leading powers of z in the two asymptotic solutions ( |20| ) near the regular-singular point z = 0. [36] 
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<fti r = 2M\s\(<f)i r — const • (f) r ) 

= -2M\s\ (p r [7 |s | \nz + z~ lsl (% + ^z + ■ ■ •) + const] . (32) 

It is straightforward to expand this expression about z = 0. This expansion yields 

i ( \ _ $ A~ s (l + afz + a+z 2 + ■•■)+&/ <t>t ln * (for s > 0), 
^ rlrj \ (1 + ar^ + a^^ 2 + •••)+ ^0" In ^ (for s < 0), 1 j 

where af and /? s j = — 2M|s|7| s are constants. 

The above analysis, based on Eq. fl28|) , explains the origin of the logarithmic term and 
determines its exact form. The calculation of 7y (and hence of f3 s {) for general I, s is 
tedious, however. It is easier to derive the explicit expression for (3 s i directly from the 
exact expression ( P4"| ) for <pi r . The series expansion of the hypergeometric function around 
(1 + z)~ l = 1 (corresponding to z = 0) may be obtained from a generating function through 
the formula 



(/!) 2 (/- |s|)!(Z + |s|)!dy' 



dy 



(34) 



(cf. Eq. (4) in Sec. 2.2.2 of |3JJ). Note that in our case y = 1/(1 + z), so ln(l - y) = 
Inz — ln(l + z). The logarithmic term in Eq. ( p3|) (which comes from the first of the above 
two ln terms) is only obtained when none of the 21 derivative operators d/dy in Eq. (|34] ) acts 
on ln(l — y). Thus, for the sake of calculating the logarithmic coefficient, we can replace the 
second factor in the right-hand side of Eq. (|34]) by 



d 



\nz = (-l)'(O!^ [(1 - vY^y' 1 - 1 ] lnz. (35) 



dy 1 

Evaluating this expression at y — 1, to the leading order in z = 1 — y, we obtain 



(H)! 



|s| lnz. (36) 



Substituting this in Eq. (|34|) , and recalling Eq. (|25|) , we obtain the desired expression for 
the logarithmic coefficient: 

^ - „ H ( :;;K)i (^)-'- + °- (37) 

In summary, we have constructed a basis of solutions to the static field equation. One of 
the basis solutions (0 r ) is simply a polynomial in 2, but the other (4>i r ) contains a logarithmic 
term. This logarithmic term will play an important role in the analysis below. Note also 
that for both s > and s < 0, (p r is smaller than 0j r in the leading order by a factor A' S L 
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V. REGULARITY AT THE EH 



By general considerations, we expect physical perturbations to be regular and smooth 
at the EH. The function \l/ s represents a perturbation in the Maxwell field tensor F a p for 
s = ±1 and in the Weyl tensor C a ^ for s = ±2. When this perturbation is expressed 
in Kruskal coordinates (|4]) (or in any other coordinates which are regular at the EH), all 
components of these Maxwell or Weyl tensors must take a perfectly regular form at the EH. 

To discuss the regularity of at the EH, it is useful to define 

& = A S V S , (38) 

and correspondingly, 

^slm = A s ^slm and psl m = A s F slm_ ^ 

It is straightforward to show, via Eqs. (0) and @, that for any s, ^ s directly represents 
a physical perturbation field which must be regular at the EH: For s = 2, ^ s is a linear 
combination of Weyl components CyaVb, where the indices a, b represent the two angular 
coordinates 8, if. For s = —2, \l/ s is a linear combination of Weyl components CuaUb- 
Similarly, for s = 1 and s = — 1 \l/ s is a linear combination of Maxwell components Fy a and 
F a u, respectively. Therefore, a necessary condition for regularity at the EH is that \l/ s be 
regular (i.e. finite and smooth). Since the spin- weighted spherical harmonics are smooth, 
%l) slm must be smooth too. 

We point out that the regularity of \l/ s is also dictated by mathematical considerations, 
as follows: If one transforms the master equation @ from to and from the original 
coordinates to the Kruskal coordinates (f|), the field equation becomes perfectly regular at 
the EH (whereas with the original dependent variable the equation is singular at the EH, 
even in Kruskal coordinates). Therefore, from the hyperbolic nature of the field equation, if 
the initial data for are regular (which we assume), no irregularity may evolve at the EH. 

Consider next the regularity of the static solutions. We assume that for any s and any 
I, m, there exists (at least) one static solution which is physically regular at the EH. For, 
if there is an external static source of a multipole I, m (and no incoming waves from past 
null infinity), the field outside the BH will be static; and we do expect this static field to 
be regular at the EH. The presence of two independent regular static solutions (for a given 
s, I, m) at the EH would violate the no-hair principle, because then all static solutions would 
be regular at the EH, including the one which is regular at infinity. We shall now show, 
however, that for any s, one of the static solutions (the solution <p ir ) is irregular. 

For s < 0, the irregularity of (p ir is obvious, because the corresponding field 4> ir = A s (j) ir 
diverges like A s . For s > the field 0j r is finite (= A ) at the EH. Yet, the logarithmic 
term implies that the solution is not smooth: The derivative of order \s\ with respect to r 
(which itself is a regular coordinate) diverges. On the other hand, for both s > and s < 0, 
the field <fr r = A s <fr r is a polynomial in z [which is proportional to (r — 2M) S for s > and 
to (r — 2M)° for s < 0], so it is perfectly smooth. 

We conclude that for both positive and negative s, <pi r is physically irregular, whereas 4> r 
is physically regular. 
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VI. LATE-TIME BEHAVIOR 



The demand for regularity of i/j slm at the EH has immediate implications to the late-time 
expansion (|Hf). Since r and v form a regular coordinate system for the Schwarzschild back- 
ground (the so-called "ingoing Eddington coordinates"), ip slm must be a perfectly smooth 
function of r and v at the EH (recall that r and v are related to U and V by an invertible 
analytic transformation). Therefore, for any k (and any s), F(! lm (r) = A s i^ /m (r) must be a 
smooth function of r: 

F s k lm (r) e C°°{R) for all k. (40) 

In section |i| we have shown that must be a static solution. The regularity of 

Fjl l J^(r) then implies that F% l J% must coincide (up to some multiplicative constant) with the 
regular static solution r . Hence, to the leading order in 1/v, we obtain 

^j slm (t, r) = c <f> r (r) v- 21 - 3 + 0(*r 2/ ~ 4 ), (41) 

where Cq is constant. For the description of the late-time behavior along worldlines of fixed 
r > 2M, it is useful to re- write this expression in terms of powers of l/£:0 

^ /m (r, t) = c Mr) t~ 21 - 3 + 0(r 2l ~ 4 ). (42) 

From this point on we discuss the cases s < and s > separately. 



A. The case s < 

In this case, Fj^ is proportional to <f)~. We shall denote the proportionality constant 
by Cq, that is, 

FS = c o -0 r - ( a <0). (43) 

(Recall that the parameter k$ in Eq. ( fL3|) is so defined such that the term F^^ does not vanish 
identically. Therefore, by definition, the constant Cq in non-zero.) Note that <p~ = A' 5 ' near 
the EH.' 

Consider next the contribution from the terms k > 0. From Eq. ( [10]) it is obvious that 
for s < and for all > 0, F£ lm must be a regular function of r, which vanishes at least like 
A' s ' at the EH (like for k — 0). Hence, at late time the terms k > are negligible compared 
to the term k = 0, due to their higher negative powers of 1/v. Therefore, Eq. ([H|), which 
now reads 

r lm (r,t)^ Co -<P;(r)v- 21 - 3 (s<0), (44) 

provides a useful description of the late-time behavior not only at r > 2M but also at the 
EH. To the leading order in A, the asymptotic behavior at the EH is 

^ m (r, t) = c A' 3 v- 21 -' 3 (s < 0). (45) 



3 It should be stressed here that in this paper we assume the power index 21 + 3 derived by Price 
H for the tail at fixed r > 2M. The new information in Eq. (42) [or in Eq. (|4l|)1 concerns the 
explicit form of the radial function multiplying the inverse-power factor. 
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B. The case s > 



In this case, F£ l J£ (which must coincide with the regular static solution) is proportional 
to 4>ti i- e -> -F/Sd = c o 4>t w ^ n c o 7^ 0) hence F£ l J% oc A near the EH. However, for k > 
(for which F^ lm is not the static solution), the only obvious constraint on the functions F^ lm 
is the regularity of F£ lm (r), Eq. (PPf). This regularity criterion allows the terms k > to be 
proportional to A _s (and, as we show below, at least the term k — 1 is indeed proportional 
to A~ s ). Due to this A~ s factor, at the EH the 0(v~ 2l ~ i ) term in Eq. ( |4"1~D dominates the 
0(v -2 ' -3 ) term, which is only proportional to A . Therefore, for s > 0, Eq. ( |4~1~D does not 
provide a useful description of the asymptotic behavior at the EH as it does for r > 2M. 
A correct description of the late-time behavior there must include both terms k = and 
k = 1: 

^™(r, t) = c+ 0+(r) v' 21 - 3 + F° l ™{r) v' 2l ~ 4 + 0(^- 5 ) (s > 0). (46) 
To the leading order in A, the asymptotic behavior at the EH is 

^j slm (r, t) = c+ v- 21 - 3 + c+A- s v- 21 -* (s > 0), (47) 



where is the coefficient of A~ s in Fj^. Note that Eqs. ( f4if ) and (f42"f ) still provide a 
correct and useful description of the late-time behavior along any line of fixed r > 2M. 

It is important to verify that the coefficient cf in Eq. ([|7]) is non-vanishing. This coeffi- 
cient is to be obtained from the function F^l(r) in Eq. (^). F^[ satisfies the inhomoge- 
neous equation (|14D , subject to the regularity condition (flop . The general inhomogeneous 
solution takes the form 

F^(r) = a+0+(r) + 6+0+ (r) + ^(r), (48) 

where and are constants and is a specific inhomogeneous solution. Using the 
Wronskian function W(r) given in Eq. pTj), we can express 0^ as 



Mr) = Mr) J dr - <j> ir {r) J dr. (49) 

For s > it is convenient to re-express this inhomogeneous solution in the form 
6- h (r) - [ r dr> ['dr" #Wr + fr") W(r>) S^r") 

Mr) -J dr U mr')] 2 W(r»)A(r>r m 



which is easily obtained from Eq. fl4"9"| ) by first substituting for ir , using Eq. (|2"8"|), and 
then integrating the resulting expression by parts. The form ([5(]) is advantageous as it only 
involves the homogeneous solution which has a simple polynomial formfj The source 



4 In Eq. (|50|) we have not specified the lower limit of the integration over r'. Changing the value 
of this limit amounts to adding a regular solution oc <^T, which is equivalent to re-defining the 
coefficient af in Eq. (|48|). Note, however, that the choice r' = 2M as the lower integration limit is 
forbidden, as the integral is not defined in this case. 
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term S\ is to be calculated from Eq. (|16|) with Ffe_i(r) = F (r) = Cp 0+(r). This yields, to 
the leading order in A, 

S x (r) S 16M 3 sk 4A-\ (51) 



In view of Eqs. (|23|) , (P7|), and (0), we find that the integrand in Eq. (p0|) is given, to the 
leading order in A, by 16M 3 s/c c^~A ,s-1 (r")/A s+1 (r / ). Performing the double integration, we 
obtain (to the leading order in A)[] 

Mr)=4Mk c+lnz + O(A°). (52) 

By substitution in the inhomogeneous equation (0), one easily verifies that the term In z 
in (f>ih must be multiplied by a homogeneous solution (otherwise the homogeneous operator 
D sl , acting on the logarithmic part of 0^, would yield a term proportional to In z - which 
cannot be balanced by the logarithmic-free source term); and from Eq. ([52]) and Eqs. (p3||26D 
it follows that this homogeneous solution must be proportional to Therefore, 

<j) ih {r) 4Mk c+ 0+ Inz + O(A ) (53) 



(in which the O(A ) term is logarithmic-free). We now substitute this expression in Eq. 
using the asymptotic forms (|23|) and (|33|), and keeping only the leading order (proportional 
to A~ s ) of the non-logarithmic part: 

Fjgi(r) = # [A- s + /? sI 0+ In z] + 4MA; c+ 0+ In z. (54) 

Note that the coefficient in Eq. ( fi8|) (which, in principle, is to be obtained by matching 
the solution to the late-time field at null infinity [jTT]] ) does not enter Eq. (|5"3]), as 0+ includes 
neither A~ s terms nor logarithmic terms. 

Now, Fi(r) must satisfy the regularity condition fl4"0|), so it cannot contain a logarithmic 
term. This dictates the value of the constant b\: 

6+ = -4Mk c+(3; l 1 + 0. (55) 

One can now identify the non- vanishing coefficient bf with the above leading order coefficient 
c~l of F\(r) at the EH. We conclude that the coefficient in Eq. (f47|) is non- vanishing. As 
a consequence, we find that on the EH itself the perturbation is dominated by the second 
term in Eq. (|4"?D , and hence it decays there like v~ 2l ~ 4 . 

In this section we have obtained the asymptotic form ([IT]) (and proved that the coefficient 
cf is non- vanishing) by a direct analysis of ijj slm in the case s > 0. There is yet another way 
to obtain Eq. (|47D, which, being somewhat outside the main course of this paper, we describe 
in detail in App. [A[ It is well known that each single one of the perturbation fields s = 1 and 
s = — 1 determines the full electromagnetic perturbation, i.e. the full Maxwell tensor F a p 



5 Note that no In 2 z, In 3 z . . . terms arise from the integration in Eq. (^): For s > the integrand 
is actually a rational function of r", analytic at r" = 2M. Hence, the integration over r" cannot 
produce a lnz' term. A term oc Inz arises only from the subsequent integration over r'. 
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(up to a trivial addition of the static Coulomb solution). Similarly, each of the perturbation 
fields s = 2 and s = —2 determines the full gravitational perturbation, i.e. the perturbation 
in the Weyl tensor (up to gauge, and up to a trivial addition of the static multipoles with 
I = and I = 1). In particular, ip 2 determines (p , and ^ 4 determines and vice versa. 
We use this fact in App. [A], where we obtain the asymptotic behavior for s > from that of 
s < (as we showed above, the latter is relatively simple, because for s < the term k = 
in the late-time expansion completely describes the late-time behavior at the EH). To that 
end, we shall use in App. [A] the well known Teukolsky-Starobinsky identities. 

VII. SINGLE FOURIER MODES 

Consider a solution to the field equation (0), having the form 

il>(r,t)=if> u (r)e- 4b *. (56) 

(The indices s,l,m, which, in fact, characterize both functions ip and are omitted here 
and below for brevity.) Each Fourier mode T/> w (r) then satisfies an ordinary equation which 
may be written as 

^ {rA s ' 2 ^) + [uj 2 + iusR{r) + V ls {r)\ (rA s ^) = 0, (57) 

in which R(r) and V/ S (r) are certain radial functions. The asymptotic form of this equation 
near the horizon is 

^ {A*' 2 ^) = (s/AM + tuf (A*/ 2 ^) . (58) 

The two asymptotic solutions at the EH are 

V£ = A°e lujr * and V>£ - A^e"^* (59) 

(where use has been made of the asymptotic relation e r * //< - 4M - ) oc A 1 / 2 ). At the limit uj — > 0, 
these two asymptotic solutions approach the two asymptotic static solutions, Eq. ( PD|) - just 
as one would expect. We shall now show, however, that for s > the role of regular and 
irregular solutions is interchanged as the limit uj —>■ is approached. 

In the case u^O, too, we expect one of the two solutions to be regular and the other 
one to be singular (for the same reasons as in the static case). We now substitute ip a and 
ip b in Eq. (|56|) , and construct the corresponding physical fields ip = A s ip (which, as was 
discussed in section 0, should be regular functions at the EH). We denote the functions ij) 
obtained from ijj® and ipti by 4> a an d respectively, and find 

i) a = A*V£e -<w * = A s e~ iu)U } i) h = A a V*e -iw * = e~ iu)V . (60) 

Recall that v is regular at the EH, but u is not (as the EH is a surface of finite v but infinite 
u). This implies that tp b is regular, but ip a is irregular. (For s < 0, tp a diverges at the EH. 
For s > 0, i\) a is finite, but its s-order derivative with respect to U is indeterminate at the 
EH, and higher-order derivatives diverge there.) We conclude that for uj ^ (and for both 
s > and s < 0), tp^ is regular and i/j* is singular. (This is a well known result; see p7|.) 

Let us now compare this situation to the static case, Eqs. ( ^Bp and (ffiS|). For s < 0, the 
classification into regular and irregular solutions is preserved at the limit u — > 0. However, 
for s > 0, the regular and irregular solutions switch role in this static limit! 
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VIII. SCALAR-FIELD TOY-MODEL 



To better understand the exchange of regular and singular solutions at the limit to — ► 
(for s > 0), it is instructive to consider a simple scalar-field toy-model. Let $ be a minimally- 
coupled, massless, Klein-Gordon test field on the Schwarzschild background. We make here 
the assumption that, in an appropriate gauge, the late-time behavior of the electromagnetic 
four-potential A a and of the linear metric perturbation h a p is qualitatively the same as 
that of a scalar field (this assumption is somewhat vague, especially because of the gauge 
ambiguity. Note, however, that at least for the behavior of metric perturbations along 
r = const > 2M lines, this assumption is verified in [|J.) Correspondingly, we would expect 
that the components F a y of the Maxwell tensor - which are made of terms like A a y - will 
qualitatively behave at the EH like For the same reason, we would expect F a u to 

behave at the EH like Recalling the way \l/ s=±1 is constructed from F a p by projection 
on the tetrad (|j) [see Eq. (§)], one intuitively expects that \l/ s=±1 will qualitatively behave 
as follows: 

^ s=1 oc A _1 $^ = ^ s= - 1 oc $ M = &=-\ (61) 

Similarly, for the case \s\ =2, one expects [in view of Eq. (|7|)] that 

^ s=2 oc A- 2 $. vv = ^ s=2 , ^ s= - 2 oc $ ;uu = ^ s= " 2 . (62) 

(For brevity, we shall focus the following discussion on the case \s\ = 1. Similar arguments 
apply to \s\ = 2 as well.)[] 

For any mode l,m of $, the two static solutions take the asymptotic forms 

$ r ^ 1 + 0(A) and $ ir = r* (63) 

near the EH (cf. J3|). Clearly, $ r is the regular mode, while $j r is singular (as — > -co at 
the EH). Let us denote the functions \E' S which correspond to the regular and singular modes 
by and ^f r , respectively. Recalling that in the static case d v = —d u = (1/2) d/dr* = 
[A/(2r 2 )] d/dr, we find for |s| = 1:0 

~ s ( A , s = +1 Ts f A" 1 , s = +1 , nA . 

*'r«{ A l] 8 = - 1 > *Sr«{ A 0,' s=-l ■ (64) 

Consider next a single Fourier mode (of a given l,m,u) $ = $ w (r) e~ luJt . The two 
asymptotic solutions of the radial function at the EH are obtained by substituting s = in 
Eq. 



6 A more sophisticated toy model would be obtained by replacing F a p or C a ^s in the definitions 
of the Newman-Penrose fields by $ ;Q/ g or $ ;Ql g 7 5, respectively. Here we adopt a simpler toy-model 
which is easier to calculate. 

7 In deriving the asymptotic form for ^>f, it has been assumed that d<fr r /dr does not vanish at the 
EH. This assumption is justified, as it is known |Ll] , |38[ ] that the EH-regular static scalar field $ r is 
nothing but Pi[{r — M)/M], the Legendre polynomial of order I (up to a multiplicative constant). 
At the EH we then have d<& r /dr oc dPi/dr = 1(1 + l)/2, which does not vanish (except for I = 0). 
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e iujr *[l + 0(A)}, $^e~ iW *[l + 0(A)]. (65) 

These two radial functions correspond to the field configurations 

$« = e "^t ^ e - iLJU [l + 0(A)], $ b = e - *"* = e~ iuv [l + 0(A)]. (66) 

Since u diverges at the EH (but v is regular), it is obvious that $ b is the regular solution, 
while $ a is singular. We shall denote the functions which correspond to the regular and 
singular modes in Eq. ( |66"D by \& a and respectively. Using Eq. (0), we find for \s\ = 1: 

~ f [1 + 0(A)], s = +1 - 6 f A- 1 e-*»[l + 0(A)], « = +1 

K \ e~ iwu [l + 0(A)], s = -1 ' \ Ae-^[1 + 0(A)], s = -1 ' 1 J 

(It is assumed here that the terms 0(A) in Eq. (|66l ) are non- vanishing, and, moreover, that 
their derivatives with respect to r do not vanish at the EH.) The construction of ^ a and \P 6 
ensures that the t-dependence of both functions is simply e~ tujt . Let us denote the radial 
parts of these two functions by ^%{r) and ^ (r), respectively; that is, 

vj> a (r, t) = e~ iwt ^(r), &(r, t) = e^tf* (r). (68) 

For both cases s — 1 and s — — 1, one thus finds 

^ oce^ r *[l + 0(A)], #*oc A- s e-^ r *[l + 0(A)]. (69) 



A comparison of Eq. fl64j) to Eqs. (|23|j26| ), and of Eq. (p9|) to Eq. (|59j) , reveals that for 
both cases = and w ^ 0, and for both s = 1 and s = — 1, the actual asymptotic 
form of both the regular and singular solutions agree with that obtained from the scalar- 
field toy- model. In particular, in the case s = —1, at the limit uo — > the regular uj ^ 
solution vf^ approaches the regular static solution \^* =_1 (and the singular w ^ solution 
^ approaches the singular static static solution \^* T =_1 ), whereas in the case s = +1 the 
regular and singular solutions interchange at the limit u — > 0. 

Our toy model provides a simple intuitive explanation for the difference in the role of the 
regular and singular solutions in the static and uj ^ cases. The key point is the relation 
between the two basis solutions of the scalar field itself, i.e. Eqs. ( |S3|) and ( |65| , |66D . In the 
static case, there is a "small solution" $ r and a "large solution" <3>j r . Naturally, the "small 
solution" is the regular one, and the "large solution" is singular. On the other hand, in 
the case w ^ both radial solutions in Eq. fl6"5|) are of the same magnitude. In this case, 
the fundamental difference between the two basis solutions is that, at the leading order, 
one of them (<3> a ) depends solely on u, and the other one ($ b ) depends solely on v. We 
can therefore refer to the two radial functions $ a , and $ b , as the "w-solution" and the u v- 
solution", respectively. Since v is regular at the EH and u diverges, the "v -solution" <3>^ is 
regular and the "w-solution" is singular. 

Now, the functions (which presumably represent the functions \l/ s ) are obtained in our 
toy model by differentiating $ with respect to u or v (depending on the sign of s) . Consider 
first the u ^ case (in which the two basis solutions are classified as a "f-solution" and 
a "it-solution"). When the operator d v acts on $, it naturally yields a large outcome for 
the "^-solution" , and a small outcome for the "w-solution" . On the other hand, when the 
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operator d u is applied, it yields a small outcome for the u v -solution", and a large outcome 
for the "w-solution" . Since the "^-solution" is regular and the "w-solution" is singular, we 
arrive at the following conclusion: For \[> s=_1 (which is associated with the regular 

solution is the smaller of the two basic solutions. However, for ty s=1 (which is associated 
with <&, v ), the regular solution is the larger of the two basic solutions. 

On the other hand, in the static case we have a "large solution" and a "small solution" 
(instead of a "f-solution" and a ll u solution"). The differentiation of the "large solution" 
with respect to either u or v yields a function \& s which is larger than that obtained from 
the differentiation of the "small solution" . Therefore, in the static case, for both s > and 
s < the regular solution is the smaller of the two basis solutions. 

The interchange of the regular and singular s > solutions in the transition from w ^ 
to u = may still look somewhat mysterious, because the limit u — > is a perfectly regular 
limit of the differential equation (57). The mystery may again be resolved with the aid of 



our scalar-field toy model. Let us re-write the regular uj ^ solution for $ [Eq. fl5E|)l in a 
somewhat more explicit form, 

$ fc ^e-^[l + c(w)A + 0(A 2 )]. (70) 

We assume that c(uj) is continuous and non-vanishing at the limit uj — > 0. We now calculate 
\I/^ =1 from this regular solution, via Eq. fl6l]) , keeping the leading-order in A separately for 
terms proportional to to and for terms proportional to uj°: 

* a=1 = e-^l-iuA- 1 + (4M)~ 1 cH] [1 + 0(A)]. (71) 
Restricting attention to the limit uj — > and to the leading order in A, we obtain 

^ s =i ^ e -^(_ zcjA -i + Co ) ; (72) 

where Co = (4M) _1 lim c(u). Equation (|7^) explains the change in the asymptotic form of 

the regular s — 1 solution from A -1 in the case uj ^ to A in the case u — 0. 

On the other hand, when the same calculation is carried out for s = —1, one obtains 
from Eq. (|6T|) 

r^^-coAe"™ (73) 

for the regular solution (for small uj). Thus, the regular solution for s = —1 is proportional 
to A for both u> ^ and u = 0. 

Our toy model also allows us to obtain the late-time behavior for both s > and s < 
directly from that of the scalar field. The late-time behavior for (a mode I, m of) $ near the 
EH is known to be J7,n 



^ <5>\(r) v~ 2l - ?J , (74) 

where the radial function $J,(r) is the regular static solution for the mode l,m. Equation 
(|61|) now yields at the EH 

^=-1 ^ _ A $i v -2l-3 ( 75 ) 
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and 



where <3> = & r (r = 2M), and $ x = (8M 2 ) -1 (d<I>J./dr) r= 2M- Compare these results to Eqs. 



(@D and (| 

So far we have implemented the toy model for the case \s\ — 1 only. The calculations 
in the case \s\ = 2 are straightforward too, though they are somewhat more tedious. We 
shall merely point out here that all the expressions we have derived for \|/ s=±1 are extendible 
to # s=±2 , and may be used to explain the various features of ^ s ±2 — e.g. the asymptotic 
behavior of the regular and singular solutions for both w ^ and uo = 0, and the late-time 
behavior at the EH. It should be emphasized that the late-time power index of ty s=2 at the 
EH is 2/+4 (and not 2/+5, which might naively be anticipated due to the two ^-derivatives in 
the definition of this function). The reason is that, the second-order covariant differentiation 
in Eq. ( p2| ) involves the differentiation of the afline connection. The easiest way to evaluate 
\|> s=2 is via the Kruskal coordinates (which at the EH minimize the connection's effect). One 
then finds that $ ; y oc v ~ 2l ~ A /V, and the next differentiation with respect to V then yields, 
at the leading order, $ ; yy oc v~ 2l ~ 4 /V 2 , i.e. & ;vv oc v _2 ' -4 . 

Finally, we point out that Eq. ([72]), which was derived within the framework of the scalar- 
field toy model, may also be derived for the realistic field \P S=1 , if Eq. (|73|) is assumed, using 
the Teukolsky-Starobinsky identities (the application of which is described in the Appendix). 
More explicitly, let us write the asymptotic behavior of the monochromatic s = — 1 field at 
the EH, to the leading order in A, as 

^=" 1 = o^jAe'™, (77) 

and assume that a(u) is non-vanishing at the limit u — > 0. Then, applying the Teukolsky- 
Starobinsky identities, one can easily obtain for the corresponding s — +1 field 

# s=1 oc e-^i-tuA' 1 + const ■ A ) (78) 

(for small oj). 



IX. A KERR BLACK HOLE 

The above analysis of the Schwarzschild case has immediate implications to rotating 



black holes as well. In a forthcoming paper the late time expansion will systematically 



be applied to the Kerr case, in order to determine the late-time behavior of external pertur- 
bations. Here, we shall use the above methods and considerations to derive the decay rate 
of s fields along the Kerr EH (many of the details are left to Ref. ||23|). 



In the Kerr case, the Master equation is fully separable only in the frequency domain, 
by writing 

* usim (i, r, 9, <p) = S^iO) e^e-^V^r), (79) 

where (t, r, 8, (p) are the Boyer-Lindquist coordinates, and S^ m (8) e %mLp are the spin-weighted 
spheroidal harmonics |[37||. The behavior of the radial function ip wslm (r) is then governed by 
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the well known Teukolsky equation [37|. The two asymptotic solutions of this equation at 
the EH are 



where 

A = r 2 - 2Mr + a 2 , (81) 

M and a are, respectively, the mass and specific angular momentum of the black hole, r* is 
defined by dr^/dr = (r 2 + a 2 ) /A, and 

^+ = (82) 
+ 2Mr+' v ; 

with r + = M + (M 2 — a 2 ) 1 / 2 being the r-value of the EH. [Compare the asymptotic solutions 
( jSOf) and (|59f) in the case a = 0.] We shall consider only a BH background with < \a\ < M 
(the extremal | = M, requires a separate treatment). 

In the Kerr case, the coordinate ip goes singular at the EH. Transforming to the regular- 
ized azimuthal coordinate 

(p + = (p — fl + t (83) 

(see Sec. 58 in fTTH), and substituting the solutions fl80|) in Eq. (|79|), we obtain the field 
configurations associated with the two asymptotic solutions: 



#^ m (t, r, 0, <^+) = S^ m (0) e im ^+ A- s e~ i(w - mf7+) ^, (84b) 

where 

= S s J m (9) e im *e- iuJt ^ m {r). (85) 

It is straightforward to extend the regularity criterion of Sec. [V] to the Kerr case: Here, 
too, one finds that at the EH the variable 

& = a s ^ s (86) 

must be a perfectly smooth function of the (regularized) coordinates (exactly for the same 
reasons described in Sec. |V| for the Schwarzschild case; see also f37j). 



For the application of the late-time expansion we must verify which of the above two 
asymptotic solutions is physically regular at the EH. Teukolsky |37j asserted that the regular 



solution is ^j,. This is obvious from the oscillatory dependence of \l/ a on u (and of on 
v) — as we have discussed in the Schwarzschild (w ^ 0) case. One must recall, however, 
that this simple classification breaks down whenever w — mQ + = (in which case the above 
oscillatory factors in u and v degenerate to 1). In this case the classification is more involved. 
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We shall now restrict attention to the static^] case w = (the case w ^ is not required 
for the analysis below). In this case, the asymptotic solutions ( fg4a| , l§4b|) become 

#£=M"»( rj e, <p + ) ^ Y slrn {6, fi + )A°e imn+u , (87a) 



qu=o,.im^ ^ ^ y s '™(# ; <^ + ) A-e*^", (87b) 

where y s/m denotes the spin-weighted spherical harmonics. Teukolsky's assertion concerning 
the regularity of ^ is now valid for m ^ only, and we still need to find out what is the 
regular asymptotic behavior at the EH for m — 0. 

Fortunately, at this point we can directly apply the results from the above analysis of 
the Schwarzschild case, for the following reason. Let us define 

_ r-r + 
2VM 2 - a 2 

The relation between A and z is A = 4(M 2 — a 2 )z(z + 1) (note that as a — > both A 
and 2! coincide with their above Schwarzschild definitions). One can now verify that for 
m = 0, to = the master equation takes exactly the form of Eq. ( F2T|) p3[ . Therefore, the 
two static solutions in the Kerr case are exactly <p r and <pi r defined above (viewed as functions 
of z). We already know that the solution <p r (like A s r ) is a perfectly regular polynomial of 
z (and hence of r), whereas the solution <j) ir (like A s ir ) includes a term proportional to lnz 
and is hence irregular at the EH. 

Let us summarize the above results concerning the regularity of static (i.e. w = 0) modes: 

(i) The case m ^ 0: For both s < and s > 0, the regular solution is ipb (just as in the 
w 7^ Schwarzschild case). The field associated with this regular asymptotic solution is 
given in Eq. ( |87b| ). 

(ii) The case m = 0: For both s < and s > 0, the regular solution is r - just as in the 
static Schwarzschild case. The field associated with this regular solution is 

Y s > l > m =°(e)<f) r {z), (89) 

with the function 4> r (z) given in Eq. (^), and its asymptotic behavior (for both positive 
and negative s) given in Eq. (H). [Note that in terms of the limit m = of Eqs. (|87a|J87b|) , 
Eq. (|8~9"D conforms with \l/ a for s > and with ^ for s < 0.] 



The above results (whose detailed derivation is given in J23|) are summarized in Table 
I. This table displays the asymptotic form of the regular and irregular static modes for the 
various possible values of s,m. 

After we have discussed the regularity features of the static solutions, we are in a position 
to analyze the decay rate of the late time tails along the EH, using the late-time expansion. 



8 Throughout this section, which deals with a Kerr background, we refer to the i-independent 
solutions as " static" , in a slight abuse of the usual terminology. (We prefer to use here the term 
"static" instead of "stationary" in order to simplify the terminology and preserve the semantic 
analogy with the Schwarzschild case). 
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As we mentioned above, the master equation for the Kerr background is only separable in 
the frequency domain. Since the late-time expansion is carried out in the time domain, we 
cannot take advantage of the full separability of the field equation. The dependence on ip 
is still separable via e l7nip , however, and without loss of generality we shall consider a field 
\l/ sm of a single m (the overall perturbation field will be obtained by a superposition of all 
m values). To deal with the dependence on 6, we proceed as follows: We first perform 
the late-time expansion of the full perturbation field \l/ sm , and then, for each term k in this 
expansion, we separate the angular dependence by decomposing into spin-weighted spherical 
harmonics. The full decomposition thus takes the form 



•qj sm = ^ 



k=0 



slm i 

k \ 



,l=lo 



v- ka ~ k , (90) 



where Iq is the minimal value of I allowed for the mode m, s in question, that is, Iq = 
max(\m\, \s\). The parameter ko is defined here to be the dominant late-time power index 
of \fr sm along lines of constant r > r + ; Namely, it is determined by the multipole I which 
has the slowest decay at r = const > r + . Note that by this definition, k is independent of I 
(unlike the Schwarzschild case, in which the late-time expansion was implemented for each 



mode I, m in separate). An investigation of the late-time decay at fixed r indicates 
that generically the dominant multipole is the one with the smallest I allowed, i.e. I = Iq, 
and its decay rate (at fixed r > r+) is v with all other multipoles decaying faster. This 

means that generically ko = 21q + 3, and also, the term k = includes only one multipole, 
I = l (that is, Fj£?o vanishes for all I > l ). 



When the expression ( |9"0|) is substituted in the master equation |!7|], one finds that the 
radial functions F£ lm still admit equations of the form (ini). However, in the Kerr case the 
source term Sf} m involves also contributions from other values of I. (Actually, the source 
term Sf m couples a multipole I to multipoles I ± 1,/ ± 2.) Still, one finds that, as in the 
Schwarzschild case, S^ m depends only on functions F k i with k' < k [p^j . In particular, the 
function Ff^ has no source term, so it satisfies a closed homogeneous equation, which is 
just the static field equation. This structure allows one to solve for all unknowns Fjf m in 
an inductive manner, starting with the functions F^q. (The situation here is analogous to 
that of a scalar field in a Kerr spacetime, analyzed in Ref. P2"|.) 

As we have just explained, the function Fk =0 must be a static solution of the master 
equation. Furthermore, the regularity arguments discussed in Sec. (V| (for the Schwarzschild 
case) are applicable to the Kerr case as well, and imply that Fk =0 must be the regular static 
solution. The decay rate of the late-time tail along the Kerr EH now follows immediately 
from the above discussion of the regular static solutions, as we now describe. 

For m 7^ 0, the regular static solution is ?/v Since i/jb( r ) has the maximal amplitude 
allowed by the regularity criterion, the terms k > 1 (being proportional to v~ k "~ h ) will be 
negligible. Therefore, the late-time tail at the EH will be proportional to ipb(r)v ~ h °, and for 
both s < and s > we shall have (to the leading order in A and 1/v) 

$ sm oc Y s > l °> m (e, <p+) e imn+v A' 3 [v~ ko + O(v- ko ~% for m ^ 0. (91) 

(The oscillatory factor e tmn+v has already been observed in the scalar- field case 

Note that the angular dependence in this expression, as well as in Eqs. (0,^) below, only 

includes the multipole Iq: As was mentioned above, generically F£ l J^ vanishes for all I > l . 
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For m = 0, the situation is just as in the Schwarzschild case: In the case s < 0, the regular 
static solution r is proportional to A~ s . Since this is the maximal magnitude allowed by 
the regularity criterion, the terms k > 1 will be negligible in this case too (just as in the 
case m/fl above). The late-time tail at the EH will therefore be proportional to <p r (r)v~ ko , 
and we obtain 

^ sm ocY s ' l0 ' m=0 (9)A- s [v- k0 + O(v- k °- 1 )}, form = 0,s<0. (92) 

However, in the case s > (yet with m — 0), the k = term is proportional to the regular 
static solution <p r oc A , whereas the k — 1 term diverges like A~ s . Therefore, at the EH 
the k = 1 term dominates, and the late-time asymptotic behavior near the EH is given by 

m sm oc Y s ' l °' m=0 (9) [v~ k ° + cA~ s y-* " 1 + 0(v- k °- 2 )}, for m = 0, s > 0. (93) 



In Ref. |y] we verify (by calculating Fk=\) that the coefficient c is non- vanishing, and also 
that the term proportional to A~ s v^ k °^ 1 includes the multipole I = Iq only. Note that for 
m = 0, s > 0, the decay rate along the EH is v whereas in all other cases it is v~ k °. 
(The decay rate along lines of constant r > r + is v ~ k ° in all cases.) 

As was mentioned above, generically ko = 21q + 3. Thus, the most dominant m-modes 
are those with \m\ < \s\. For these m-modes Iq = \s\, so fco — 2|s| + 3. Therefore, in the 
overall perturbation field made of the superposition of all m-modes, the decay rate along 
the EH is generically ti~ 2 l s l -3 . Recall, however, that for s > fields the axially-symmetric 
component (m = 0) decays faster along the EH, like v - 2 N -4 . 

Note that the discussion is this section only deals with non-extremal, \a\ < M, Kerr 
BHs: the extremal case needs to be considered separately. (Note, for example, that Eq. 
is not valid in the case a = M.) 



X. CONCLUSIONS 

To summarize our results, it is most useful to refer to the physical variables = A s ^ s 
defined above. These variables are natural, because, as we mentioned above, at the EH they 
are proportional to the regular Maxwell components F a y or F a u, or to the regular Weyl 
components C a vbv or C a ijbu (for s = 1, s = —1, s = 2, and s = —2, respectively; here, a, b 
stand for the two regular angular coordinates, i.e. 9, <p in the Schwarzschild case and 9, tp+ 
in the Kerr case). 

For the Schwarzschild case, we find from Eqs. ( pE5| ) and ( fT7D that along the EH, each 
mode l,m of the physical field ^ s decays at late time with the leading-order tail 



^P slm (v) = const • v - 2l - 3+p , for s < 

^slm/ v \ ^ congt . for s > o 



Schwarzschild BH, (94) 



where the constants are generically non-vanishing. Here, p = 1 if there exists an initial 
static multipole I, or p = in the absence of such a multipole. (In the calculations above we 
assumed p = 0, but the extension to the p—\ case is trivial.) For the comparison with the 
Kerr case below (in which a full angular separability of the late-time tails is not possible), 
it is useful to re-write Eq. (|94"D in terms of the field \l/ sm , which is the part of \l/ s including 



22 



all multipoles I for a given m. ^ sm is dominated by the minimal multipole I allowed for the 
values s, m, i.e. Iq = max(\m\, \s\). Equation ( |94] ) yields 



qsm ^ congt . y s ,i ,m^ ^ v -2i -3+p^ for s < Q 
* sra ^ const • K sA) > m (0, ^) u -2« -4+p ) for s > 



Schwarzschild BH. (95) 



This holds for both axially-symmetric and nonaxially-symmetric modes. 

For the (nonextremal) Kerr case, we find from Eqs. fl9~T| ) — (^) that the physical field 
^,sm £ Qr a gj ven value of m decays along the EH according to 

S const • y s ^> m (#, £ + ) e imQ+,; w - 2 '«- 3 +f, Kerr BH, m ^ 0, (96) 

for nonaxially-symmetric modes, and 

q sm ^ const .Y s ' lo ' m=0 (e)v- 2l °- 3+p , fors<0] „ ^ TT 

; n ,J „ /, , r" Kerr BH, m = 0, (97) 

^ sm ^ const • F s ' /o ' m=0 (^)w- 2Z °- 4+p , fors>0j ' ' V ; 

for axially-symmetric modes. 

The overall late-time behavior of the field is dominated by the m-values which yield 
the minimal possible value of Iq, i.e. the values — \s\ < m < \s\, for which Iq = \s\. Thus, 
the overall decay rate of \l/ s is obtained from Eqs. (^,^,^) by substituting l — > \s\. Note, 
however, that the angular dependence of the overall late-time field (and, in the Kerr case, 
also the oscillations in v) will be obtained by a superposition over all m values in the range 
— \s\ < m < \s\. 

The late-time behavior of an s > field along the EH of the Kerr black hole displays 
two important differences between the axially-symmetric and nonaxially-symmetric modes. 
First, the modes m^O oscillate along the horizon's generators according to e imn+v , whereas 
the mode m = decays monotonically. Second, the modes m ^ (with \m\ < s) decay 
like v~ 2s ~ 3 , whereas the mode m = decays along the EH like v _2s_4 . (The first difference 
applies also for s < fields, but the second one is special to s > 0.) These two differences 
lead to an interesting consequence: The overall s > late-time field oscillates along the 
generators of the Kerr EH (because the non-oscillatory mode m = decays faster than the 
oscillatory modes m/0). On the other hand, for s < fields the overall late-time tail at 
the Kerr EH is a superposition of oscillatory m^O tails and a monotonic m = tail (which 
decays at the same rate). 

The above difference between the s < and the s > fields in the power-law indices 
of the tails along the EH has never been reported before (as far as we know), even in the 



Schwarzschild case. In Ref. |26[ Hod attempts to calculate these tails in the Kerr case, but 
his analysis yields no difference between the s < and s > fields, even for m = (see, 
however, the note added). In what follows we briefly explain what seems to be the reason 
for the incorrect result in Ref. ||26| . 

Hod uses the correct asymptotic radial solution for uj ^ modes near the EH, which for 
m = reads ijJu{r) = C{oj)/Sr s e~"' lwr * (see Eq. (39) in Kl fT). He then continues by assuming 



9 Note that in j2(| a different dependent field variable is used, which is A s / 2 times the Newman- 
Penrose fields that we use in the present paper. 
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that C(u) is cu-independent. This assumption, however, is invalid in the case am = 0,s > 0, 
where C(u) oc to. In the Schwarzschild case (for any m), this result was demonstrated in 
the framework of our toy-model in Sec. IVTl (] [see Eq. ([72])], and can also be verified for the 
actual s > fields by using the Teukolsky-Starobinsky identities [see Eqs. ( |77| , [78D 1. In the 
Kerr case (for m = 0), the situation seems to be the same.Q If the correct form, C(u) oc oo, 
is used in the case m = 0, s > 0, then, Eq. (40) in |26| correctly yields a tail smaller by 1/v 
than the s < tail. 

The asymptotic behavior of the various physical fields along the EH is important for 
understanding the dynamics of these fields inside the black hole: One can naturally view 
the field value at the EH as initial data for the black-hole's interior. (This is most naturally 
implemented within the framework of the characteristic initial- value formulation.) Of special 
importance is the case of gravitational perturbations (s = ±2) of the Kerr background. In 
this case, evolving the perturbation from the EH to the future, determines the gravitational 
field (and hence the spacetime geometry) inside the black hole, up to the inner horizon. The 
infinite blue-shift of the gravitational perturbations leads to a curvature singularity at the 
inner horizon ||. Knowing the late-time behavior of the perturbation along the Kerr EH 
enables one to analyze in detail the structure of this singularity . 

Note added. — After this paper was submitted, Hod presented the correct result for s > 
fields (the one derived in Sec. [IX| above), in a recent manuscript |28 . 



APPENDIX A: DERIVATION OF ^ s>0 FROM * s<0 USING THE 
TEUKOLSKY-STAROBINSKY IDENTITIES 

In this appendix we obtain the asymptotic form of \l/ s>0 at the EH of a Schwarzschild 
BH from that of \1> S<0 , using the Starobinsky-Teukolsky identities. By this we shall recover 
Eq. (f47|), and show that the leading-order coefficient does not vanish. (This was already 
verified in Sec. flVTl ) by a direct calculation of ty s>0 from the field equation for s > 0.) 

The Teukolsky-Starobinsky identities [EO] relate the perturbation fields ty s>0 and \l/ s<0 . 



In the case of the Schwarzschild background, these identities take the form 

(a s - iub s d t ) V£(r) = Vl s [^- S (r)] (s > 0), (Al) 

where ip^{r) is the radial Fourier mode introduced in Eq. (p6|), a s is a non- vanishing constant, 
b s is a constant which vanishes for s = 1 but not for s = 2, and V is a differential operator 
given by 

V (r) =d r - ioor 2 A-\ (A2) 



10 As was shown in Sec. IX, for m = 0, s > (in Kerr) the regular solution switches from ipj, to 
tp a as u vanishes (just like for s > in the Schwarzschild case). This switching seems to indicate 
the vanishing of C{uj) as to — > 0, as discussed in Sec. VIII . [Note that for m / 0, s > 0, the regular 
solution is V'b for both u> = and w^O (we are only interested here in small uj values near uj = 0, 
so we can now assume u — m£l+ ^ 0), indicating that in this case C(w) is non- vanishing as uj — > 0. 
The same holds for s < (and any m).] 
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(There also is an analogous identity for transforming from s > to s < 0.) If we now apply 
the inverse Fourier transform to Eq. (Al), we obtain for the time domain function ip slm (r, t) 



[the one introduced in Eq. fllTD] 

(a s + b s d t ) r lm = V s ^- S ' lm } (s > 0). (A3) 

where T> s is the differential operator 

V s (r, t) = (d r + ^A^dt) 23 = (2r 2 A- 1 9„) 2s . (A4) 

Here, the f-derivative is taken with fixed u, and the t-derivative is taken with fixed r. 

Before making use of this identity to study the late-time tails, let us briefly discuss its 
application to the static solutions. In the static case, Eq. ( |A3|) reduces to 

a s ^ slm = {d r ) 2s ip~ s ' lm (s > 0, static). (A5) 

Consider first the application of this identity to the regular static solution <p~ (namely, we 
take s > and if)~ s > lm = <p~). Since <p~ is a polynomial in r, the right-hand side must 
be a polynomial too. Since the outcome must be a static solution for s > 0, it must be 
the polynomial static solution, i.e. 0+ (up to some constant). This confirms our previous 
conclusion, namely, that for s >0 the regular static solution is 0+ and not <pf r (i.e. the one 
proportional to A and not to A _s ). 

Next, consider the application of the identity ( |A5| ) to the other, irregular, static solution 
0~ . When the differential operator (d r ) 2s acts on the logarithmic-free terms in the right-hand 
side of Eq. (|3"3"D (2nd row), it yields a regular polynomial, as before. However, when applied 
to the logarithmic term, it yields two types of terms: (i) logarithmic terms, proportional to 
A - these are obtained if the derivative operator never acts on the factor ln(z). (ii) Non- 
logarithmic terms proportional to negative (as well as positive) powers of A: These terms 
are obtained when one of the operators d r acts on the factor ln(z). The most dominant 
negative power is obtained when d r acts on ln(^) on its \s\ + 1 operation, which yields a 
contribution proportional to A~ s (all other contribution are of less negative powers of A). 
One can identify the terms (i) and (ii) with the second and first terms, respectively, in the 1st 
row on the right-hand side of Eq. (|3~3| ) . Note the crucial role played by the logarithmic term 
in 4>~ (despite the fact that it only appears in a sub-dominant term proportional to A' s l): 
Without this logarithmic term, the operation in Eq. (A"5) would have yielded a perfectly 
smooth function of r, proportional to A . 

We shall now apply the (time-domain) Teukolsky-starobinsky identity ( |A3|) to the s < 
late-time field (H). For the consistency of the notation, we shall assume that s > [as 
dictated by the notation of Eq. (|A3|) 1, and therefore re- write Eq. fl44|) as 



*Jj- s > lm (t, r) = Co #T(r) v' 21 ^ (s > 0). (A6) 



Applying the differential operator T> s to the right-hand side of Eq. (|Aq), we obtain three 
types of terms: 

(i) a term in which the derivative operator d v never acts on v~ 2l ~ 3 , 

(ii) a term in which the derivative operator d v acts on v~ 2l ~ 3 once, 

(iii) terms in which the derivative operator d v acts on v~ 2l ~ 3 more than once. 
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Consider first the term (i). As was discussed above, the operator t> s transforms <p r into 
<f>r- Hence the term (i) is nothing but 

Cq v- 21 - 3 V s {(j)~) = const • v~ 2l ~ 3 • 0+ (A7) 

(with a non- vanishing constant). 

The term (ii) is proportional to v~ 2l ~ A . We must keep this term, however, because its 
radial function will appear to diverges at the EH [whereas the radial function in the term 
(i) is regular]. The regular static solution 0~ is, at the leading order, proportional to A s 
(recall that now s > 0). The most divergent contribution is obtained when the differential 
operator d v acts on v~ 2l ~ 3 on its \s\ + 1 application (with the other contributions smaller by 
factors of A). This contribution yields 

d v (v~ 2l ~ 3 ) (9 r ) s - 1 [2r 2 A- 1 (9 r ) s (A s )] = 

[2(2/ + 3)(-l) s s!(s - l)!(2M) 2s+1 ] A~ s iT 2 ^ 4 [1 + 0(A)]. (A8) 

(Recall that when acted on a purely radial function, 2r 2 A~ 1 c\, = r 2 A _1 <9 r * = d r .) 

The terms (iii) decay as v~ 2l ~ 5 or faster. The radial functions involved in these terms do 
not diverge faster than A~ s (in fact, they diverge even slower). Therefore, these terms may 
be neglected. We find that 

V s {*fj- S > lm ) = (i) + (ii) = (co v- 21 - 3 + c 1 A- s tr 2 '- 4 ) [1 + 0(A) + 0(l/v)}, (A9) 

where both constants c"o, C\ are nonvanishing and proportional to Cq . 

So far we considered the contribution to T> s (ip~ s,lm ) from the term of ip~ s ' lm proportional 
to v~ 2l ~ 3 , which is given in Eq. ( |A6| ) . This leading-order term corresponds to the term k = 
in Eq. (|13|) . We now consider the contribution of the k > terms. From Eq. (|40|) we learn 
that for s < all functions F£ lm are smooth functions of r, which vanish at least like A' s ' 
at the EH. One can now easily analyze the contribution of each term k > in the same 
way the dominant k = contribution was analyzed above. Again one obtains contributions 
which are analogous to the terms (i),(ii) or (iii) below, except that these contributions are 
now multiplied by an extra factor v ~ k . The contributions from all k > terms of ip~ s > lm can 
therefore be neglected, and we are left with Eq. (A9). 



Once t> s (ip~ s ' lm ) is known, we can calculate ip slm via Eq. (|A3|). For s = 1, the coefficient 
b s vanishes, so ip slm = T> s (ip~ s ' lrn ) / 'a s . For s = 2, the left-hand side of Eq. ( |A3| ) includes 
the derivative operator c\. To extract ifj slm in this case, we apply the differential operator 
a s + b s d t to the right-hand side of Eq. ([T3| ) (recalling d t — > d v ), and solve for F£ lm term by 
term by matching the powers of 1/v to Eq. (|A9|) . For k = 0, this matching yields 

a s F^ = d + O(A). (A10) 

For k — 1 we obtain 

a s F° l ™ - (21 + 3)b s F° l ™ = c 1 A~'[l + 0(A)], (All) 
which, in view of Eq. ( |A10|) , we simply write as 

a s F% 1 ™ = cxA- s [l + 0(A)]. (A12) 



26 



We thus obtain 

^ slm = ( c + v - 21 - 3 + C +A- S v- 2l ~ A ) [1 + 0(A) + 0(l/v)] (s > 0), (A13) 
with = c 0) i/a s . 

Thus, relying on the late-time behavior (|44]) for s < 0, and using the Teukolsky- 
Starobinsky identities, we have recovered Eq. p7j ) for s > — with non-vanishing coef- 
ficients Cq and cl . 
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TABLES 





irregular static solution 


regular static solution 


am = 0, s > 




A u 


am = 0, s < 


A 


A~ s 


am ^ 0, s < and s > 




/\— s p irr&l+v 



TABLE I. The asymptotic forms of the physically-regular and physically- irregular static solutions at 
the EH in the Kerr case. Presented are the results for the axially-symmetric (m = 0) modes of the fields, as 
well as for its nonaxially-symmetric (m ^ 0) modes. [The Schwarzschild case (a = 0) can be read from this 
table by referring only to the results in the first two lines (which then apply to all to).] 
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